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Abstract 


The stability of the nonconservative double pendulum model of 
Ziegler is studied with emphasis placed on the effect of changes in the model 
parameters and on the presence of retardation cf the angular position of the 
follower load. Both the cases of external and internal damping are considered. 
As in previous studies, the damping is assumed to be small. 


t 


Introduction 


The stability of nonconservative elastic systems has been the sub- 
ject of much recent interest. To demonstrate some of the properties of such 
systems several authors have investigated the behavior of a simple two-di- 
mensional model consisting of an inverted double pendulum subjected to a fol- 

2 

lower type of load (Figure 1) Ziegler [1] discovered that the presence of 
linear viscous internal damping may tend to destabilize the system by lower- 
ing the value of the critical follower load. Herrmann and Jong [2,5] studied 
the behavior of this model in the case of small internal damping for a fol- 
lower load and a partial follower load. 

The results of these studies indicate that damping and the type of 
loading may have significant effects on the stability of the system and that 
considerable care must be exercised in the modeling of the damping mechanism; 
Herrmann and Jong have recommended further study of other damping models and 


loadings. 


A type of loading which might occur in some systems and which is 


probably present in most experimental studies is that of a retarded load, 

that is, a load that acts in a prescribed manner in relation to the system, 

but which is delayed by the inherent time lag of the mechanism which produces 

it. Kiusalaas and Davis m investigated the stability of the double pendulum 

* 

model subjected to a retarded follower force with a constant time lag for the 
case of no damping. 

The purpose of this study is to examine the stability of the stand- 
ard double pendulum model when subjected to a retarded follower load with small 




time lag but with the inclusion of slight damping, both internal and external. 
It is found that the presence of damping greatly alters the critical load of 
the system, an expected result. The effect of the time delay on the follower 
load is shown to be significant and, in most cases, destabilizing. Yet, for 
a few cases the time delay raises the critical load above that for the stand- 
ard follower load, an unexpected result which is significant for the inter- 
pretation of experimental data. 

A second aspect of this study is the determination of the effects 
of changes in the model parameters. Previous works have rather consistently 
studied the double pendulum with mass ratio m j/ m 2 * 2, length ratio 

2 ■= 1 and spring constant ratio e^/cg * 1 (Figure 1). It is of in- 
terest, for design and analysis purposes, to have a knowledge of the effects 
of these ratios on the stability of the system. 

/ 

The Equations of Motion ^ 

" " 0 

Consider the double pendulum of Figure 1, consisting of two rigid, 
weightless bars of lengths and j}g which carry concentrated masses m^ 

and mg, respectively. The generalized coordinates <t>^(t) and <t>g(t), with 
t the time, specify the angles between the vertical and the two bars, and are 

assumed to remain sufficiently small. The linear elastic restoring moments 

... 

• , « 

c^4>^ and C 2' < * > 2” < * > 1^ anc ^ the linear viscous internal damping moments b^^ 
and ^2^2"‘^1^ are assume< 3 to act at the two hinges, respectively, with b^ 
and bg small but not both zero. 

The load P is applied at the free end at an angle 0(t) with 
respect to the vertical, and it is assumed that 


f 
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0(t) * ♦g(t-T), (1) 


with t §: 0, a constant time lagj hence, the load acts at time t in such a 
manner as to be tangential to the position of the upper bar of the pendulum 
at time t-T. The special case t = 0 reduces P to the standard follower 
load. It is assumed that t is small. 

Consideration of Lagrange* s equations in the form 


d ,&S . A bD 


dT t dV 

o<t>i 


i -1,2 


(2) 


with the kinetic energy, dissipation function, potential energy and generalized 
forces given by 


T " 2 m l^l + 2 m 2^1^1 + ^2^2^ 9 


D - \ Vr + i VV 1 /- 


v - i V? + s ° 2 <V*i> 2 ’ 


( 3 ) 


\ * F M V e) ' 

Q 2 = 


immediately yields the two linearized equations of motion 


I 


(m 1 +m 2 ) J e^ 1 + (b^bg)^ + + ^jf^2 + 


-bg4> 2 - Cg<>g + P^© * 0, 


VlVx ’■ Vl ‘ c 2 *l + W 2 + b 2*2 + 


+ (Cg-P^g)^ + P^g© * 0. 


Given the linearity of the equations, assume a solution of the form 


♦ 4 (t) ■ V ' 1 * l. 2 


and recall that, from Equation (1), 


6(t) = A^e 0 *. 


Since the time delay t is assumed to be small, let 


e(t) A 0 (l-CDT)e a>C . 


Substitution of Equations ( 5 ) and ( 7 ) into ( 4 ) immediately yields the condi- 
tion for existence of a nontrivial solution as 


|i ( 8+l)n 2 +(B 1 +B 2 )n+€+l-nF (l0 2 -BgQ-l<+lF( 1-ofl) 


« 0 , ( 8 ) 


wfl —BgJi - 1 


Cl +Bgf2+1 -Raft 




I 


where the parameters of the system have been appropriately nondimensionalized 
and are defined by 


f - i 0 - “ ■ 7 - 


& 2 m 2 


b i , , „ . i . . m i c i 

B i 373 > 1 * 1 > 2 < 8 ' ;r > « = “ • 

1 i 2 (c 2 m 2 ) 1 / 2 *2 m 2 C 2 

Expansion of the determinant (8) yields the frequency equation in the form 


p 0 ^ + + P 2 fi2 + P 5 fi + P4 s °t 


( 10 ) 


from which, through the familiar Routh-Hurwitz criteria, the conditions for 
asymptotic stability are obtained as 


p^ > 0, 1 = 0, » • • , 4 

x S ( Pl P 2 - P 0 P 3 )P 3 - P4P1 > O' 


( 11 ) 


The Standard Model 


Consider now the case for which 


H ■ 1, 6 = 2, € * 1 


m 


which is the "standard” model discussed in References [1-4]. For this particu- 
lar case the coefficients of the frequency equation take the form 
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P r 2 
0 

P 1 = B 1 + 6B 2 “ 2aF 

P 2 = 7 - 2F + B^Bg - CkB-jF - 2aBgF (15) 

= B x + Bg + OF(F-5) 


Under the assumption of slight damping, the case of physical interest, and of 

small time delays the three last terms in Pg are neglected since they are 

• • 
higher order quantities. By assuming that Bg £ 0 and defining the ratios 


P = 


5 

B, 


V 


v a 2 

X “ 1” * TT » 
B 2 *2 


(14) 


the stability criteria (11) yield, for this standard model, the following four 
inequalities on the load parameter F; ' 


P< M> 

f < Qr 9 
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F<£ , 


F(VF) < 


P+1 


2y^f^ - r(i2+2p+i4r)F 5 + r(T8+i5P+20r)F 2 + 


-[ 2p 2 +l4p+12+104y+19rp]F + 4p 2 +53P+4 > 0. 


(15) 


These inequalities determine the critical load! the vertical position of the 
double pendulum is an asymptotically stable equilibrium state for loads below 
the critical load and an unstable equilibrium state for higher loads . 
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The values of these critical loads were obtained with the use of a 
digital computer and are plotted in Figures 2a and 2b. In Figure 2a the criti 
cal load is plotted versus the ratio of damping coefficients p for various 
fixed values of y, the measure of the time delay. The case y = 0 corres- 
ponds to the result obtained in [2]. It is clear that for high values of y 
an increase in p always leads to an increase in the critical load; for low 
values of y, however, the critical load increases up to a certain value of 3 
and then decreases slightly thereafter: the optimum choice of p depends on 

the magnitude of y (for y = 0, the maximum critical load F * 2.086 occurs 
at P s 11.071 [2]) , It is also noted that for some values of y a sudden up 

ward jump in the critical load takes place as p is increased beyond a cer- 

■ * ’ 

tain value. Mathematically, this occurs when two of the real roots of the 
quartic (15) become complex conjugates. 

, The critical load is plotted versus the time delay parameter y 
in Figure 2b for various fixed values of the damping ratio 'p. The presence 
of a time delay usually produces a destabilizing effect on the system, with a 
larger time delay producing a lower critical load. Since physically the im- 
portant problem is that with small damping, it is clear that the delays will 
have to be small for any degree of stability; hence the reason for the assump- 
tion on the smallness of t. Note, however, the surprising result that in the 
case p = 5 a small time delay stabilizes the system; at y ** 2, for example, 
the critical load is 4 per cent higher than for the standard follower load 
(r = o). 
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E ffect of External Damping 

Consider now the standard model of Figure 1 but with the internal 
damping of the previous section replaced by viscous external damping which 
acts on the masses m^ and m^ (but not on the bars) with coefficients b^ 
and bg, respectively. The dissipation function in (5) becomes 


d . I v* + 1 WV : 


(16) 


Following the same procedure as in the previous section, the fre* 

i ■ 

quency equation (10) is obtained with coefficients 


P. 


B 1 + 2B 2 “ 2aF ’ 

1 - 2F + B^Bg - ©R^F, 

B x + 5B 2 - 2BgF + OF(F-5), 


= 1# 


from which the stability criteria 


( 17 ) 


F < 


P+2 

2 r 




F(2+3r-YF) < P+5, 

2r 2 F^ - r(4+2p+i4r)F 5 + (30r+i5rP+20r 2 +4p)F 2 + 
-(2p 2 +20p+8+44y+19yp)F + 4p 2 + 25P + 16 >0, 


(18) 
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are obtained. The critical loads determined from these equations are depicted 
in Figure 3a, as plotted against p, and in Figure 5b, as plotted versus y. 

It is to be noted that in the case of external damping the presence 
■f a time delay is always destabilizing: an increase in y causes a decrease 
in the critical load (for fixed p). The behavior of the pendulum under 
changes in the parameter p again depends on y. For y > 1/2 the critical 
load increases with increasing p, but for y < l/2 the critical load first 
increases and for values beyond a certain value of (3 then decreases. 

For a standard follower load (y s 0), the maximum critical load 
occurs for (3 = 2, where it attains the value F = 2.0 86 (the same maximum 
value as for the case of internal damping). This value happens to be the 
value F [2] of the critical load of the model for zero damping. Hence, 
external damping on the double pendulum exhibits a "destabilizing effect" 

similar to that caused by internal damping: the addition of slight linear 

0 4 - 

/ 

viscous external damping lowers the critical load below the value F (ex- 

, 6 

x 3 

cept for p * 2, when the critical load equals F ). 

Effect of Parameter Changes 

Thus far, only the "standard" model with 5*2, € = 1 and p = 1 
has been considered. It is desirable to investigate the effect of the mass 
ratio 5, the spring constant ratio e and the length ratio p on the sta- 
bility of the model with internal damping. 

The procedure used for this purpose is identical to that used in 
the previous sections. The results of this analysis are shown in Figures 
4-6. Figures 4a,b depict the critical load for various values of 5 with 

% ~ ------- ; _ : i r • - — r -, r -Tr-- r- r - . - 

It must be noted that this "destabilizing effect" of external damping does not 
occur for a continuous cantilevered column subjected to a follower load [5]. 


€ = 1 and n s 1; the results are shown as a function of p in Figure 4a 
with a constant time delay parameter T = 1, and as a function of y in Figure 
4b for a fixed damping ratio p * 1. It is noted that the critical load in- 
creases with increasing mass ratio 6. 

The effect of various spring constant ratios e is shown in 
Figures 5a, b. It is observed that for p > 1 the critical load decreases 
with increasing e, but that the opposite effect occurs for small values of 
■ p. Figures 6a, b depict the effect on the critical load of the length ratio 
p. for fixed values 5 = 2 and € = 1. 

From these graphs it is noted that this model is extremely sensitive 
to parameter variations, and that therefore considerable care should be exer- 
cised in inferring conclusions about continuous elastic systems ,xom results 
obtained from a simple two-degree-of“ freedom model. 

Concluding Remarks 

" ■ "■ ■" ■ 

It has been shown that the presence of load retardation, a phen- 
omenon almost inevitable in experimental studies, has a significant effect 
on the stability of the Ziegler model, the effect being in general desta- 
bilizing. The importance of the nature of the damping mechanism, already 
noted by other authors, and the high sensitivity to parameter variations 
of the model have been studied and general quantitative results have, been 
obtained, 

Acknowled gment 

The authors wish to thank Dr. C. M. Strauss for carrying out the 


computer analysis. 


11 

References 

[1] H. ZIEGLER, Die Stabilitatskriterien der Elastomechanik, Ingenieur- 
Archiv., vol. 20, 1952, pp. 49-56. 

[2] G. HERRMANN and I. C. JONG, On the Destabilizing Effect of Damping in 
N onconservative Elastic Systems, Journal of Applied Mechanics, vol. 52, 
Trans. ASME, vol. 87 , Series E, 1965 , pp. 592-597. 

[5] G. HERRMANN and I. C. JONG, On Nonconservative Stability Problems of 

Elastic Systems with Slight Damping, Journal of Applied Mechanics, vol. 
55, Trans. ASME, vol. 88 , Series E, 1966 , pp. 125-155. 

[4] J. KIUSALASS and H. E. DAVIS, On the Stability of Elastic Systems under 
Retarded Follower Forces, to appear in Int. J. of Solids and Structures, 

[5] R.H.PLAUT and E. F. INFANTE, The Effect of External Damping on the Sta - 

bility of Beck* s Column, to appear in Int. J. of Solids and Structures. 



Captions for Figures 


Fig. 1 Double pendulum model 

Fig. 2a, b Standard model, internal damping 

Fig. 5a, b Standard model, external damping 

Fig. 4a,b Effect of changes in mass ratio 

Fig. 5a, b Effect of changes in spring constant ratio 

Fig. 6a, b Effect of changes in length ratio 


















